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ABSTRACT 



Context. L1448-mm is the prototype of a low-mass Class protostar driving a high- velocity jet. Given its bright H 2 spectra observed 
with ISO, L1448-mm is an ideal laboratory to observe heavy water (HDO) emission. 

Aims. Our aim is to image the HDO emission in the protostar surroundings, the possible occurrence of HDO emission also investi- 
gating off L1448-mm, towards the molecular outflow. 

Methods. We carried out observations of L1448-mm in the HDO(lio-ln) line at 80.6 GHz, an excellent tracer of HDO column 
density, with the IRAM Plateau de Bure Interferometer. 

Results. We image for the first time HDO emission around L1448-mm. The HDO structure reveals a main clump at velocities close 
to the ambient one towards the the continuum peak that is caused by the dust heated by the protostar. In addition, the HDO map 
shows tentative weaker emission at = 2000 AU from the protostar towards the south, which is possibly associated with the walls of 
the outflow cavity opened by the protostellar wind. 

Conclusions. Using an LVG code, modelling the density and temperature profile of the hot-corino, and adopting a gas temperature 
of 100 K and a density of 1.5 x 10 s cm -3 , we derive a beam diluted HDO column density of ~ 7 x 10 13 cirr 2 , corresponding to a 
HDO abundance of ~ 4 x 10~ 7 . In addition, the present map supports the scenario where HDO can be efficiently produced in shocked 
regions and not uniquely in hot corinos heated by the newly born star. 

Key words. ISM: individual objects: L1448 — ISM: molecules — stars: formation 



1. Introduction 

Water can be considered as the most important molecule in the 
oxygen chemistry (after CO), influencing in particular the over- 
all chemical composition of the gas involved in the process of 
star formation. In addition H2O is the most abundant compo- 
nent in the grain mantles and plays a fundamental role in the en- 
l— 1 ergy balance of star-forming regions (e.g., Ceccarelli et al. 1 19961 
O ! Kauf man & Neufeld [19961 Bergin et al. [19981 van Dishoeck & 
r 7^ ' Blake [19981 . In starless regions water abundance is very low (< 
J> , 10" s , e.g., Bergin & Snell 2002 ), while it increases to the abun- 



dance of CO (< 10~ 4 ) near protostars (e.g. Liseau et al. 119961 
Harwit et al. [1998). Several mechanisms have been proposed 
for H2O formation (see the recent work on H2 18 by J0rgensen 
& van Dishoeck 120 101 and reference therein): (i) heating from 
forming stars, (ii) jet-driven shocks, or (iii) shocks in accretion 
disks. In particular, water emission becomes even more impor- 
tant when studied in low-mass protostars, which are possibly as- 
sociated with solar-type protoplanetary systems. 

Unfortunately, the H2O low-energy lines are almost invis- 
ible from the ground due to very strong atmospheric absorp- 
tion, and only under very special conditions are they observed 
towards the brightest sources (Cernicharo et al. 1990). As an ex- 
ample, H2O emission has been detected towards the low-mass 
protostar NGC1333-IRAS4B by different facilities, namely ISO, 
Odin, SWAS, Spitzer and, very recently, Herschel (Maret et al. 
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2002; Bjerkeli et a l. 120091 Bergin et al. 2003; Franklin et al. 
[19891 Watson et al. 120071 Kristensen et al. l2TJT0b . These obser- 
vations, with their relatively large apertures, revealed water from 
different components, namely the protostar, the outflow and the 
extended ambient medium. Recently J0rgensen & van Dishoeck 
(120101) imaged IRAS4B in the mm-range with the PdB interfer- 
ometer. They detected H2O emission associated with the pro- 
tostar, probably coming from the inner 25 AU of the accretion 
disk, which calls for further observations of low-mass objects. 

Given the difficulty of detecting water, several alternatives 
have been pursued. The most natural one is to observe rare 
isotopologues like HDO, which suffer much less telluric ab- 
sorption and have some lines at sub-mm and mm-frequencies. 
Furthemore, if HDO is formed on the surfaces of grains before 
the heating due to the radiation from an accreting protostar (e.g., 
Parise et al. 2005a and references therein), then its emission can 
also be used to trace the past history of water formation. Heavy 
water observations have been successfully carried out towards 
high-mass star-forming regions (Plambeck & Wright 1987; Jacq 
et al. [19901 Gensheimer et al.|T9"9"6l Helmich et al. [19961 Pardo 
et al. 1200 II Comito et al. 2003), but they are still very sparse. In 
Orion, the emission arises from warm gas (Tk > 150 K) and the 
abundance of HDO is strongly enhanced (> 10 3 ) with respect to 
ambient values, most likely because of recent HDO evaporation 
from dust grain mantles (see also Bergin et al. 2010, in the HIFI 
A&A special issue). On the othe hand, the number of studies of 
HDO towards the lower-mass star-forming regions is even less. 
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Stark et al. (12004ft and Parise et al. (12005 ai ) observed HDO emis- 
sion with the JCMT and IRAM single-dish antennas, towards the 
hot corino heated by the IRAS 16293 protostar. In this case the 
measured deuteration ratio (3 10~ 2 ) supports models where wa- 
ter is produced in the gas phase at low temperature and then stor- 
aged on the grains (e.g., Roberts et al. 2000). On the other hand, 
if H2O is formed in shocks, followed by condensation out on 
the grain mantles, the deuteration ratio is predicted to be lower 
than 10~ 3 . It is still missing a study of HDO emission through 
interferometric high angular resolution images, which is needed 
to distinguish the hot-corino emission from that associated with 
more extended outflows. 

In this Letter, we report for the first time a high-resolution 
image of HDO at 80.6 GHz towards the protostar L1448-mm. 
We show that the emission at velocities close to the systemic 
one comes from the hot corino, we also report tentative HDO 
emission offset from the protostar, which is possibly associated 
with the cavity of the outflow red lobe. 

2. Observations and source selection 

L1448-mm, located at ~ 250 pc from the Sun, can be considered 
as the prototype of a low-luminosity (~ 8 L G ) Class protostar. 
The newborn star is driving a highly collimated molecular jet 
consisting of well separated extremely high-velocity (> 50 km 
s~') molecular bullets, well detected in SiO and CO (Guilloteau 
et al. |1992al Bachiller et al. 1995 ). Lower velocity CO emission 
delineates a biconical cavity, possibly created by the passage of 
the jet, which is entraining the surrounding gas through succes- 
sive bow-shock well detected in H2 (e.g., Davis & Smith 1995; 
Dionatos et al. 120091 Neufeld et al. l20U9l . L1448-mm has been 
selected because of its bright H2O spectra that are observed with 
the ISO satellite, which lead to a large water abundance, 5 10" 4 , 
towards the protostar (Nisini et al. 120001 1. The large ISO beam 
(80") does not allow one to exactly localise the HiO-rich gas, 
but evidence is given that the high water abundance is main- 
tained on a high level along the southern red-shifted part of the 
outflow. 

The HDO observations of L1448-mm were obtained in July 
2009 with the IRAM Plateau de Bure Interferometer (PdBI) in 
France. Four tracks in the five-element D configuration (base- 
lines from 24 m up to 94 m) were used, for a total time of 
~19 hours. A region of about 60" around the protostar has 
been mapped. The HDO(lio-ln) line at 80578.295 MHz was 
observed with both 20 and 40 MHz bandwidths and a corre- 
sponding spectral resolution of ~ 0.13 and 0.25 km s , respec- 
tively. Two other units with 320 MHz bandwidth have been used 
to measure the continuum. We used the same setup for verti- 
cal and horizontal polarisation to gain in sensitivity. Amplitude 
and phase were calibrated by observing 0130-171, 0234+285 
and 0333+321, whereas the flux density scale was derived by 
observing 3C84 and 3C 454.3, with an uncertainty of ~ 20%. 
Images were produced using natural weighting, and restored 
with a clean beam of 6'.'29 x 4'.'87 (PA=123°). Note that we 
observed HDO(lio-ln) emission with a deep integration (230 
minutes on-source) using the 30-m IRAM single-dish antenna 
and a HPBW of 31" (see Sect. 3). 

3. Results: the HDO maps 

Figure[T]shows the contour map of the HDO( lio-l 11) integrated 
emission (contours) towards L1448-mm. The magenta contours 
trace the unresolved 3.7 mm continuum emission with a peak 
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Fig. 1. Contour map of the HDO(l 10-1 11) emission (black con- 
tours) integrated between -0.3 and +10.7 km ss~' superimposed 
onto the continuum image at 3.7mm (magenta). The HDO con- 
tour levels range from 3<x (33 mJy km s _1 beam -1 ) to 5<x by lcr. 
Continuum contour levels range from 30 (~ 20<x) to 90% of the 
peak value (16 mJy), and identify the driving protostar L1448- 
mm. The filled ellipse in the lower left corner shows the synthe- 
sised PdBI beam (HPBW): 6'.'29 x4'.'87. The black arrow points 
in the direction of the jet (PA=-21°) identified by Guilloteau et 
al. (1992) through SiO emission. The dashed lines identify the 
sketch of the limits of the conical blue-shifted (+6, +8 km s~') 
cavity observed in CO(1-0) at PdBI by Bachiller et al. (1995). 



flux of 16 mJy centred at a J2 ooo = 03 h 25 m 38 s .9, <5 J2 ooo = +30° 44' 
05'.' 1, in agreement with previous observations (e.g. Guilloteau 
et al. U992bl Bachiller et al. 119951 The HDO structure reveals a 
main clump (hereafter called A), (i) located close to the contin- 
uum peak due to the dust heated by the protostar and (ii) sligthly 
elongated in the S-E direction, i.e. the direction of the red-shifted 
outflow lobe. The offset seen in Fig. Q] between the the clump 
A and the continuum (< 2") is marginally significant, since the 
position uncertainty is around 1". Although we cannot exclude a 
possible association of the clump A with the outflow (see below), 
we conservatively assume that clump A is tracing the hot-corino 
around the protostar. Beside clump A, the HDO map shows ten- 
tative emission from a weaker position (here called clump B) 
located at ~ 8" (^ 2000 AU at 250 pc) from the protostar, again 
towards the southern direction. No hints of HDO emission have 
been detected towards the N-W direction, where the blue-shifted 
molecular lobe is located. 

Figure|2]shows the HDO(l 10-I n) spectra, obtained from the 
PdBI maps by integrating over the whole area where heavy wa- 
ter emission is detected (upper panel) and over clumps A and 
B (middle panels). In particular, a Gaussian fit of the clump A 
spectrum leads to a peak of 49 mK (in 7mb scale), observed with 
a S/N = 5, at +6.0 ± 0.4 km s _1 , i.e. slightly red-shifted with re- 
spect to the systemic velocity. The integrated intensity is 286 
± 50 mK km s -1 , while the FWHM linewidth is 5.7 + 1.0 km 
s . The spectra, smoothed to 1.4 km s resolution to increase 
sensitivity, do not allow us a careful analysis of the line profile. 
Bottom panel of Fig.|2]reports the HDO spectrum as observed at 
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Fig. 2. HDO(lio-ln) spectra (in T% scale) derived from the 
PdBI image and integrated over the whole emitting region (up- 
per panel), and the area of the two structures called A and B 
(middle panels). The spectral resolution is 1.4 km s . Bottom 
panel reports the HDO(lio-ln) low-spectral resolution (5 km 
s" 1 ) spectrum as observed with the 30-m IRAM antenna. The 
ambient LSR velocity (+4.7 km s" 1 ) is marked by the vertical 
dashed line. 



low-spectral resolution (5 km s" 1 ) after 230 minutes on-source 
using the 30-m IRAM antenna. The HDO line peaks at ~ 6 mK 
(1 cr = 0.2 mK). The line flux is 305 mJy km s" 1 , showing that 
the PdB image collects about 61% (187 mJy km s" 1 integrated 
over the 30-m beam) the of the flux observed with the IRAM 
single-dish. Given the calibration uncertainties of the PdBI im- 
age (~ 20%) and of the 30-m measurement (< 15%), the most 
plausible solution is that a weak extended structure has been fil- 
tered out by the interferometer. 

We can compare the present HDO image with the CO(l- 
0) map obtained at the PdBI with a spatial resolution of ~ 3" 
by Bachiller et al. (1995 see their Fig. 5). The dashed lines 
in Fig. Q] identify the position of the walls of the conical out- 
flow cavity as observed at red-shifted velocities from +6 to +10 
km s" 1 ). The black arrow shows the direction of the protostellar 
jet identified by Guilloteau et al. (1992) through SiO emission. 
Considering also the HDO spectrum of clump B, which shows 
emission at the same velocities, it is tempting to associate the 
HDO emission tentatively observed offset from the L1448-mm 
protostar for the first time with the cavity opened by the jet and 
well observed in CO. Interestingly, one possible origin for the 
H2O masers at 22 GHz observed at high-angular resolution stud- 
ies towards low-mass protostars is its association with conical 
shells (Marvel et al. 120081 1. which may be caused by propaga- 
tion of large bow-shocks that entrain ambient medium at long 



transverse distances from the jet but which could also represent 
walls of evacuated cavities produced by a wide-angle wind (see 
e.g., Shang et al. 120071 and references therein). The water maser 
spots require extremely high-density (~ 10 7 — 10 9 cm 3 , Elitzur 
et al. 1 19891 Kaufman & Neufeld[1996) conditions, which could 
be found in compressed shocked layers at the interface with the 
surrounding medium. Indeed, Cernicharo et al. (1996) reported 
a spectrum at 183 GHz observed with the 30-m IRAM antenna, 
which was owing to the H20(3 13-220) transition. The H2O spec- 
trum is composed by several narrow red-shifted components, 
resembling that classical maser spectral pattern. The brightest 
feature is at ~ +9 km s _1 , i.e. close the velocity observed with 
the present HDO emission, but the whole pattern covers a broad 
range of velocities, questioning a possible association of the 183 
GHz emission with the cavity walls. 

In conclusion, the present HDO map supports the presence of 
water in the cavity walls and suggests that HDO can be found not 
only in heated hot cores/corinos, but also in shocked material. 



4. The HDO abundance inside and outside the hot 
corino 

To derive first the HDO column density in clump A and then 
the abundance, we used an LVG code with the collisional coeffi- 
cients from Green (1989), in a plane parallel slab (where because 
the line is very thin optically the exact geometry is negligible). 
We considered a gas temperature in the 50-120 K range and a 
density between 1 X 10 7 cm" 3 and 1 x 10 9 cm" 3 . Figure 3 shows 
the beam-averaged HDO column density corresponding to the 
signal observed in the direction of clump A as a function of the 
gas temperature for various densities. Within the 50-120 K in- 
terval of gas temperatures, the HDO total column density (N lot ) 
varies very little as function of the gas temperature, and mostly 
depends on the gas density. At 100 K, varying the density from 
~ 10 8 to ~ 10 7 cm -3 , corresponds to a factor ~ 4 in the HDO col- 
umn density. At higher densities, the LTE solution is approached 
and the dependence on the gas density is not as strong. 

If one assumes that the unresolved HDO emission towards 
the source is caused by the sublimation of the water ice mantles 
when the dust temperature reaches 100 K, one can derive the 
abundance of the sublimated heavy water. We used the density 
and temperature profiles of the envelope of LI 448 -mm derived 
by J0rgensen et al. (2002) in order to evaluate the sizes and den- 
sity of the sublimation region and its H2 column density. The 
sublimation region is predicted to have a radius of ~ 25 AU, 
equivalent to a diameter of about 072. At that radius, the H2 den- 
sity is equal to 1.5 x 10 s crrT 3 , giving an H2 column density of the 
sublimation region of about 1 . 1 x 10 23 cirT 2 . The beam-averaged 
HDO column density derived from the Figure 3, assuming a den- 
sity of 1.5 X 10 s cnT 3 and a temperature of 100 K, is ~ 7 x 10 13 
cm" 2 . Thus, taking into account the dilution factor of the PdBI 
beam, the HDO abundance in the hot corino results in ~ 4x 10~ 7 . 
The main uncertainty in this value is due to the uncertainty in the 
radius, density and H2 column density of the sublimation region. 
A density a factor 3 larger would decrease the estimated HDO 
abundance by less than a factor 2, whereas a factor 3 lower den- 
sity would increase the abundance by about a factor 4. Taking 
the worst scenario, i.e. that the emission arises in a region of 575 
size, and using again the structure by J0rgensen et al. (2002), 
the emitting region would have a radius of ~ 700 AU, where the 
gas temperature is ~ 22 K and the density is ~ 1.5 x 10 6 cm" 3 . 
The HDO column density (~ 1 x 10 16 cm" 3 ) and the abundance 
(~ 4 x 10" 7 ) would be outrageously large for such a cold tem- 
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Fig. 3. Beam-averaged HDO total column density corresponding 
to the ( 1 lo—l n) signal observed in the direction of the hot corino 
(clump A) as function of a gas temperature T gas between 50 and 
120 K and for a density between 1 x 10 7 cm" 3 and 1 x 10 9 cm" 3 . 
An LVG code in a plane parallel slab has been used (see Sect. 
4). 



perature, where water would be largely frozen onto the grain 
mantles, and we accordingly exclude this possibility. 

The linewidth may bring some additional information. The 
observed HDO linewidth is around 6 km s _1 . Again, follow- 
ing J0rgensen et al. (2002 ), if the emission originates in the as- 
sumed 25 AU, the mass of the central object is about 0.1 M . 
Alternatively, if the central object has a mass of 1 M , the emis- 
sion must originate in a region at about 100 AU, where the gas 
temperature is ~ 50 K, the density is ~ 2 x 10 7 cm" 3 , and the H2 
column density is 6 x 10 22 cm" 2 . The beam-averaged HDO col- 
umn density would be 5 x 10 14 cm" 2 (Fig. 3) and consequently 
the HDO abundance ~ 5 x 10" 7 . 

In summary, the present observations suggest that the HDO 
abundance is ~ 4 x 10" 7 and likely originates in a region whose 
radius is not much larger than about 100 AU. The most probable 
hypothesis is that the HDO emission observed in clump A is 
associated with the hot corino, and the sublimation of the icy 
mantles. 

We also compared our observations with the theoretical pre- 
dictions by Parise et al. (2005b), who carried out a theoretical 
study of the HDO emission in L1448-mm. Their Fig. B2, up- 
per right panel, shows the dependence of the HDO flux as func- 
tion of the HDO abundance in the outer envelope (where the 
dust temperature is less than 100 K) and in the hot corino. The 
curves are degenerate, in the sense that the observed signal is 
obtained either when (1) the HDO abundance is ~ 2 x 10~ 8 in 
the outer envelope and lower than ~ 5 x 10~ 8 in the hot corino 
or (2) the HDO abundance is ~ 5 x 10" 7 in the hot corino and 
lower than ~ 3 x 10" 9 in the outer envelope. The present in- 
terferometric observations suggests that the emission observed 
towards L1448-mm originates in the hot corino, thus solving the 
degeneracy. Also, note that the abundance derived in our simple 
analysis agrees very well with the abundance derived using this 
more sophisticate modelling. 

As discussed in the introduction, HDO can be a useful probe 
of the corresponding main isotopologue, H2O, which is largely 
undetectable from ground. The study of the water deuteration in 
a similar solar type protostar, IRAS 16293-2422, by Parise et al. 
(2005a) suggested that the HDO/H2O ratio in the sublimated ices 



is around 3%. When comparing L1448-mm and IRAS 16293- 
2422 with other sources where the molecular deuteration of 
formaldehyde and methanol was measured, it appears that both 
L1 448-m m and IRAS 16293-2422 are typical sources (Parise et 
al. [20061 . It is then plausible that the HDO/H 2 in L1448-mm 
is similar to what was found in IRAS 16293-2422, namely ~ 3%, 
implying a H2O abundance of about ~ 1 x 10" 5 . This is some- 
what similar to the water abundance estimated in IRAS 16293- 
2422 by ISO observations (Ceccarelli et al. 2000; Crimier et al. 
2006 ), and troubling because it would imply either an abundance 
of the water ices lower than usually assumed, ~ 10" 4 , or that 
not all ices have been sublimated. New observations by the just 
launched Herschel satellite will soon confirm or disregard this 
result and clarify the situation. 

The HDO emission tentatively observed outside the hot 
corino and presumibly associated with the outflow red-shifted 
cavity observed in CO is only a factor 2 weaker than that ob- 
served towards the protostar (see HDO spectra of clumps A and 
B reported in Fig. 2). Therefore, as for the hot corino region, this 
HDO signal corresponds to a HDO column density of about 10 14 
cm -2 if the density is around 10 8 cm" 3 and the emission fills up 
the PdB beam. It is difficult to estimate the abundance of HDO 
because of the lack of available estimates of the H2 emission 
in the cavity. Even though we are unable to quantify how much 
water has been sputtered or sublimated from the grain mantles, 
clearly one or more shocks are responsible for the HDO spot that 
is coincident with the cavity wall. In other words, HDO emission 
cannot be just limb-brightned emission: water must have been 
extracted from the grain mantles by some energetic process. As 
previously noted by other authors in other contexts, models of 
bow shock naturally would explain the morphology and the pres- 
ence of shock "far away" from the high-velocity jet. 

5. Conclusions 

The present high spatial resolution map of L1448-mm allows us 
to image for the first time HDO emission at 80.6 GHz around a 
Class protostar with a spatial scale of 5'.'5, i.e. ~ 1400 AU at 
250 pc. We have shown that the bulk of the emission at veloci- 
ties close to the ambient one likely comes from the hot corino, 
whereas a tentative detection of HDO emission is present at ~ 
2000 AU from the protostar, and it is possibly associated with 
the walls of the outflow cavity previously observed in CO. These 
findings support the idea that HDO can be efficiently produced 
in shocked regions and not uniquely in hot cores heated by the 
newly born star. 

We derived the HDO abundance in the hot corino with 
the reconstruction of the density and temperature profiles by 
J0rgensen et al. d2002l l. Adopting a gas temperature of 100 K 
and a density of 1.5 x 10 8 cm" 3 , a HDO column density of ~ 7 
x 10 13 cm" 2 was derived, corresponding to a HDO abundance 
of ~ 4 x 10" 7 . If we assume that the HDO/H 2 in L1448-mm is 
similar to what was found in IRAS 16293-2422 (~ 3%), a H 2 
abundance of about 10" 5 is derived. Only future estimates of the 
H2 column density will allow us to give an estimate of the HDO 
abundance in the outflow. 
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